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Negative differential conductivity in far-from-equilibrium 
quantum spin chains 
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Abstract. - We show that, when a finite anisotropic Heisenberg spin- 1/2 chain in the gapped 
regime is driven far from equilibrium, oppositely polarized ferromagnetic domains build up at 
the edges of the chain, thus suppressing quantum spin transport. As a consequence, a negative 
differential conductivity regime arises, where increasing the driving decreases the current. The 
above results are explained in terms of magnon localization and are shown to be structurally stable 
against breaking of integrability. 



Quantum transport properties of low-dimensional ma- 
terials are currently an object of intensive theoretical and 
experimental research [1-3] . Low-dimensional systems are 
interesting for theoretical investigations, as they admit 
ordering tendencies, leading to collective quantum states 
that are difficult to realize in three-dimensional systems. 
Experimentally, unconventional transport properties have 
been reported, including unusually high thermal conduc- 
tivity in quasi-one dimensional magnetic compounds and 
ballistic spin transport in magnetic chains (see Ref. [1] and 
references therein). Understanding the transport proper- 
ties of such low-dimensional strongly correlated systems is 
a challenging open problem. So far, most of the theoret- 
ical studies concentrated on the close-to-equilibrium situ- 
ation by using the linear-response formalism [4-6], while 
almost nothing is known about the physics of such systems 
far from equilibrium. On the other hand, new quantum 
phases and interesting physical phenomena may appear in 
the far-from-equilibrium regime [7]. 

In order to drive an interacting quantum system far 
from equilibrium, it is necessary to strongly couple it to 
some macroscopic reservoirs. Theoretical description of 
such situation usually relies on a master equation approach 
for the density matrix of the system, where the dissipa- 
tive term depends on the coupling of the model to the 



reservoirs. An analytical treatment of such situations in 
non-trivial many-body systems is generally unfeasible [8], 
while numerical simulations are highly demanding and, as 
far as we know, only very few examples are discussed in 
the literature [9]. 

In this paper we propose a conceptually simple model of 
an anisotropic one-dimensional quantum Heisenberg spin- 
1/2 chain coupled to a pair of magnetization reservoirs. 
We consider a spin chain driven far from equilibrium, 
much beyond the linear-response regime: we show that 
a long-range spin ordering into ferromagnetic domains is 
induced, independently of the ferromagnetic or antiferro- 
magnetic nature of the spin-spin coupling. This coopera- 
tive, many-body quantum state hampers spin flips, thus 
strongly suppressing the current. Therefore, the spin cur- 
rent, which at small driving obeys Ohm's law, for suf- 
ficiently strong driving exhibits negative differential con- 
ductivity (NDC), namely increasing the driving decreases 
the current. 

A NDC effect has been already predicted and observed 
in nanoscopic objects, such as single molecules, short nan- 
otubes, or quantum dots weakly coupled to metallic elec- 
trodes (see, e.g., Refs. [10, 11] and references therein). 
Transport properties in such systems are usually studied 
by considering effective models of a few single-particle lev- 
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els. Conversely, we remark that in our case the NDC has 
a different origin and unveils a new phenomenon: it arises 
as an outcome of a beautiful interplay between coherent 
many-body quantum dynamics of the spin chain and in- 
coherent spin pumping. 

Our system is constituted by N interacting quantum 
spins 1/2, whose dynamics is described by an XXZ Heisen- 
berg exchange Hamiltonian: 
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x, y, z) are the Pauli matrices of the fc-th 
spin, and A = J z / J x denotes the xz anisotropy. Both ends 
of the spin chain are coupled to magnetic baths, that is, 
to magnets acting as reservoirs for the magnetization [4]. 
Within the Markovian approximation, which holds pro- 
vided the bath relaxation time scales are much shorter 
than the time scales of interest for the system's dynam- 
ics, the time evolution of the system's state density matrix 
pit) follows the Lindblad master equation [12]: 
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are four Lindblad operators, crjr = (a^ ± ia\)/2 are the 
raising and lowering operators, while [•, •] and {•, •} de- 
note the commutator and the anti-commutator, respec- 
tively Hereafter we use dimensionless units, by setting 
h = 1, J x = 1. The Lindblad operators L\, L2 act on 
the leftmost spin of the chain (fc = 1), while L3, L4 on 
the rightmost one (k — N). The dimensionless parame- 
ters T and pl (pr) play the role of the system-reservoir 
coupling strength and of a left (right) chemical potential, 
respectively. In other words, 2pl,r—1 6 [—1, 1] is the cor- 
responding bath's magnetization per spin in dimensionless 
units. 

Given the conceptual simplicity of our model, various 
experimental implementations could be envisaged. For in- 
stance, one could consider molecular spin wires [13, 14] 
with each boundary coupled to a magnetic impurity, where 
the desired populations {pl,r, 1 — Pl,r} of up/down spins 
can be pumped by means of electromagnetic fields. We 
choose a symmetric driving: pl,r = |(1 T so that 
p = pr — pl £ [0, 1] is a single parameter controlling the 
driving strength. When p is small we are in the linear- 
response regime, while in the limiting case p — 1 (corre- 
sponding to pl — 0,Pr = 1) the left (right) bath only 
induces up-down (down-up) spin flips. This last case can 
be thought of as two spin-polarized chains coupled to the 
system. For instance, one could have two oppositely po- 
larized ferromagnets for the left and the right baths. 




Fig. 1: (Colour on-line). Stationary spin current (j) as a func- 
tion of the driving strength p, with anisotropy A = 2 and 
system-baths coupling strength F = 4. We rescaled the y-axis 
with the system size N in such a way that the small- p Ohmic 
behaviour can be clearly appreciated. 



The time evolution of the master equation (2) has been 
numerically simulated by employing a Monte Carlo wave 
function approach, based on the technique of quantum tra- 
jectories [15, 16]. We assume that there exists a single 
out-of-equilibrium steady state and study the stationary 
spin current (j) , which is calculated by looking at the left 
bath and summing up all up-down flips minus all down- 
up flips, and then dividing by the simulation time. Due to 
the conservation of the total magnetization in the Hamil- 
tonian model (1), after the convergence time this current 
is precisely equal to the analogous quantity computed at 
the right bath and also equals the expectation value of 
3k = J x {a%°l+x ~ °fc°i+i) for any fe = 1, . . . , iV - 1. 

Quite surprisingly, for an anisotropy |A| > 1, the sta- 
tionary spin current exhibits a negative differential con- 
ductivity phenomenon, as shown in Fig. 1. Namely, in 
the linear-response regime of small p we find an Ohmic 
behaviour (j) oc p/N of the spin current (j); the current 
reaches a maximum at a certain value p*, and, further 
increasing p, it decreases, thus exhibiting a negative dif- 
ferential. At maximum driving strength p = 1, the current 
drops exponentially with TV. While the small-/i Ohmic be- 
haviour is consistent with recent linear-response numerical 
results [6], the NDC effect is completely unexpected. 

NDC is robust upon the variation of the number N 
of spins (Fig. 1) and of the system parameters T and A 
(Fig. 2) . By increasing T over more than four orders of 
magnitude, from T = 1CP 2 to T = 3x 10 2 , we find [see 
Fig. 2a)] that the driving strength p* at which the current 
reaches its maximum shifts from p* « 0.9 to p* ~ 0.3, thus 
considerably reducing the range of validity of the linear- 



1 Furthermore, the phenomenon has been hinted by one of us in 
a different context of kicked open quantum dynamics: Prosen T., 
unpublished notes available at arXiv: 0704 . 2252 [quant-ph]. 
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Fig. 2: Driving strength fi* at which the spin current exhibits 
a maximum (left) and current at /i* minus current at \i = 1 
(right), for a chain of N = 6 spins. Top panels are for A = 2, 
bottom panels for F = 4. 
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Fig. 3: (Colour on-line). Spin magnetization profiles (cr|) s 
versus the scaled spin index coordinate tjErr at A = 2, T = 4, 
for driving strengths /j, — 0.3 (empty symbols) and /i = 1 (filled 
symbols). 



response regime. The maximal current drop, measured 
by (i) M = M * - (j) M =i, is obtained for T s» 5 [see Fig. 2b)]. 
When r« 1 the "contact resistance" 1/r becomes large 
and the current is small for any value of fi. On the other 
hand, for T 1 we are in the quantum Zeno regime [17], 
where the coupling to the reservoirs is strong and freezes 
the system's dynamics, thus drastically reducing the spin 
current. 

In the bottom plots of Fig. 2, the driving strength fj,* 
and the current drop — are shown as a function 
of the anisotropy A. For |A| < 1 the negative differen- 
tial conductivity effect is absent and the linear-response 
regime (j) oc /i can be extended up to /x = 1, with an 
ideally conducting behaviour, so that the spin current is 
independent of the system size [5]. On the other hand, 
NDC is observed at any |A| > 1 [see Fig. 2c)]. Since spin 
transport along the chain is suppressed when |A| ^> 1, 
there exists a value of the xz anisotropy, |A| w 1.5 [see 
Fig. 2d)], at which the current drop is maximal. From the 
above analysis we infer that, for a chain of N — 6 spins, 
the optimal working point for the observation of the NDC 
phenomenon is T ss 5, |A| ss 1.5. 

Fig. 3 displays the stationary spin magnetization pro- 
files (cr^) along the chain, for the parameter values of 
Fig. 1. In the linear-response regime, we observe a con- 
stant linear gradient, with the magnetizations {crl N ) s of 
the two borderline spins close to the bath magnetizations 
( a L r) = 2/Z£,ij — 1 = =F/i. This behaviour is typical of nor- 
mal Ohmic conductors [18]. Interestingly, in the limiting 
case /i = 1 we notice the appearance of a stationary state 
characterized by two almost ferromagnetic domains, that 
are polarized as the nearest reservoir and whose relative 
width increases with the system size. These ferromag- 
netic regions are responsible for strongly inhibiting spin 



flips, and therefore for suppressing the spin current. We 
remark that we found no differences between a ferromag- 
netic (J z < 0) and an antiferromagnetic ( J z > 0) spin cou- 
pling, and in particular we observed the same stationary 
spin profiles. In the first case the formation of ferromag- 
netic domains may be somewhat intuitive, and indeed it 
can also be observed in the ground state of the autonomous 
XXZ chain exposed to two oppositely oriented static mag- 
netic fields applied at the chain edges [19,20]. On the other 
hand, the build up of ferromagnetic domains is is a priori 
not obvious for an antiferromagnetic coupling. Indeed in 
this last case, for the autonomous model (1), ferromag- 
netic domains correspond to a highly excited state; from 
the point of view of the Hamiltonian system, the net ef- 
fect of the baths is that of pumping energy into the system, 
while leading to a stationary state with very low entropy. 

As shown in Fig. 4, the time scale needed to reach the 
stationary state at fi = 1 is exponentially long. This is the 
key observation on which our intuitive explanation of the 
NDC phenomenon relies, as we shall qualitatively explain 
below. Since at small fi the current grows linearly, it is 
sufficient to show that the current is suppressed in the 
limiting case /i = 1 to conclude that, due to the continuity 
of (j)fj,, a region of negative differential conductivity exists. 
Therefore, we focus on the fi = 1 case. 

First, it is instructive to consider a situation in 
which the system is coupled to a single, fully polarized 
reservoir, [ij, = 0. Regardless of the anisotropy A, 
the stationary state is pure and ferromagnetic, namely 
III ■•• 4) (41 ■■■ \\, since the Hamiltonian (1) conserves 
the overall magnetization while at the left boundary of the 
chain only the lowering operator oc acts. As Fig. 4 
(inset) shows, at |A| > 1 also for the single-bath case the 
convergence of (cr|) to the equilibrium value {erf) = — 1 
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Fig. 4: (Colour on-line). Time approach to the equilibrium 
value of the spin magnetization of the leftmost spins, for fj, = 1, 
N = 12, A = 2, T = 4. From bottom to top: k = 1,2, 5. 
In the inset we plot the times t* at which the magnetization 
reaches a given fraction of the corresponding stationary state 
value: (cr|) = ^(a^) s , as a function of the spin index k (empty 
circles). The case with a single bath coupled to the first spin 
(k = 1) is also shown (filled squares). 

requires a time scale t* which grows exponentially with the 
spin-bath distance k (this sharply contrasts with the ide- 
ally conducting case |A| < 1, where t* only grows linearly 
with k). Moreover, if an up-polarized reservoir is added 
at the right boundary, at |A| > 1 the relaxation times t* 
for the spins closer to the left than to the right boundary 
do not change significantly. This implies that, in our two- 
baths model, the spin polarization is practically affected 
only by the nearest bath. 

The above results can be explained in terms of local- 
ization of one-magnon excitations. Given a ferromagnetic 
state |0) = HI ••• |), one-magnon excitations have the 

general form X)/cli c fcl^)i wnere W) = |0) describes 
the state with the A:-th spin flipped. If the XXZ chain 
of Eq. (1) has open boundary conditions, there is an en- 
ergy gap 2\J Z \ between states |1) and \N) (spin-flip exci- 
tations at the boundaries) and states |2), |3),..., \N — 1). 
Indeed one has (1| H s |1) = (N\ H s \N) = (N - 3) J z , and 
(2| Hs |2> = . . . = (N - 1| Hs \N-1) = (N- 5)J Z . Only 
nearest-neighbour spin-flipped states are coupled, with 
a coupling (k\Hs\k + l) = 2J X , therefore the Hamilto- 
nian (1) in the one-magnon basis {|1) , |2) , . . . , \N}} is a 
tridiagonal matrix that can be diagonalized exactly, in the 
limit of large N. One finds that there always exist at least 
N - 2 delocalized solutions. If |A| = \J Z /J X \ > 1, both 
for ferromagnetic and antiferromagnetic coupling, two pe- 
culiar eigenstates emerge: these "molecular orbitals" read 
\ip±) ~ 7f (IV'l) ± \^r}), where the states \^ L . R ) are cen- 
tered at sites 1 and N, respectively, with a localization 
length I m 1/ In |A|. The gap between the two correspond- 
ing energy levels shrinks exponentially with the system 



size, so that the coherent tunneling between the two bor- 
der sites requires a time scale that increases exponentially 
with N. In practice, this means that spin-flip excitations 
created at the border of the chain remain exponentially 
localized when |A| > 1, over a localization length £. 

Consider now the action of the two baths. As can be 
clearly seen from Fig. 4, intermediate states in the re- 
laxation to the stationary spin profile have two opposite- 
oriented ferromagnetic domains close to the baths. In or- 
der to enlarge them, spin-flip excitations should be prop- 
agated, through (J^(J^ +1 and a v k o v k+Y exchange couplings 
of Hamiltonian (1), across the ferromagnetic domains to 
the chain boundaries. Suppose, for instance, that we have 
the leftmost m spins down and the (m+ l)-th spin up and 
that this excitation propagates to the left bath; then the 
bath can flip this spin down, thus ending up with a fer- 
romagnetic domain with m + 1 leftmost spins down. The 
crucial point is that the one-magnon propagation is ex- 
ponentially localized at |A| > 1. This explains why (i) 
equilibration needs exponentially long time scales and (ii) 
only the nearest bath is felt. We stress, however, that neg- 
ative differential conductivity is observed after quite short 
time scales (t ~ 10 2 ), independently of the chain length. 
Indeed, it is sufficient to create a very small ferromagnetic 
region close to a bath, to strongly suppress current. For 
example, this spin-blockade mechanism can be clearly seen 
at [i ~ 0.9 — 0.95, where only a couple of outer spins reach 
magnetization values close to ±1, but still the current is 
far below the peak value at fi*. 

We finally mention that the presence of NDC is not re- 
lated to integrability of the Heisenberg model. Indeed, we 
performed numerical simulations on XXZ chains in pres- 
ence of a staggered magnetic field along the z direction, 
-Bstag(fc) = (— l) k B. This system exhibits a transition 
from integrability to quantum chaos [21], when increasing 
the field strength B. We found that the NDC phenomenon 
is insensitive to such a transition. 

We point out that our model might be relevant for far- 
from-equilibrium electronic transport through molecular 
wires that are embedded between two electrodes [22]. This 
is suggested by the fact that the XXZ chain can be mapped 
by the Jordan- Wigner transformation [23] into a tight- 
binding t — V model of spinless fermions moving with a 
hopping amplitude t = 2J X on a iV-site one-dimensional 
lattice, with nearest neighbour interaction of strength V = 
4J Z . Here the spin current of the XXZ model becomes a 
charge current, the Lindblad operators (3) are mapped 
into operators injecting/extracting electrons at the lattice 
boundaries, and |(1 + (erf)) into the charge density at 
site k. In particular, the ferromagnetic domains of Fig. 3 
imply a phase separation in the fermionic model with all 
the electrons frozen in half of the lattice thus inhibiting 
charge transport, if \V/t\ > 2, and, notably, irrespectively 
of whether the interaction is attractive or repulsive. 
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